JOB

"read in the data"
"indiv = individual identifier"
"eventday = day on which adverse event occurred"
"sta = day before the start of the observation period"
"end = day which is the end of observation period"
"exday = day marking the start of exposure"
FILEREAD[name = 'oxford.txt'] identifier = indiv, eventday, sta, end, exday; fgroups = no

"convert indiv to a factor using GROUPS with the option lmethod=*"
GROUPS [redefine=yes; lmethod=*] indiv

"ex1, ex2 etc.. mark either 1 day before 1st day of a risk\ 
period or the last day of a risk period for each individual.\ 
risk periods must not overlap, i.e. ex1 < ex2, ex2 < ex3 etc...\ 
(can give ex1 = ex2, ex2 = ex3 etc... to skip out a risk group)\
If there is no risk period give ex1,ex2 etc.. values greater than max(end)"
CALC ex1 = exday + 14
CALC ex2 = exday + 35

"ex: pointer containing the risk factor cut points."
POINTER [values = ex1,ex2] ex

"exc: labels for risk factors corresponding to the cut points in ex.\ 
The first and last value of exc should be the baseline level, 0" 
VARI [nvalues=3] exc
READ exc
0 1 0 :

"t: age cut points to mark the end of an age group or day before the start of\
a new age group. Do not include the beginning of the first age group or the end\
of the last age group. The minimum length of t is 1 (for 2 age groups). If\
no age groups are required choose a value greater than or equal to max(end)"
VARI [nvalues=1] t
READ t
547 :
 
"Put any covariates into a pointer named 'covariates'. All covariates are\ 
assumed to be factors (if not alter 'GROUPS [redefine=yes] cov[]' in sccs.gen)."
"POINTER [values = ] covariates"

"Delete anything that is no longer needed"
DELETE [redefine = yes] exday

"fit: what explanatory variables to fit (this should include 'age_group': the age group\ 
factor and 'exposure_group': risk factor). Covariates can be fitted as interactions with\ 
exposure_group. Covariates will be labelled cov[1], cov[2] etc.. in the order they were put into\
the 'covariates' pointer"
TEXT [VALUES = 'exposure_group + age_group'] fit

OPEN 'sccs.gen'; CHANNEL=4; FILETYPE=input
INPUT 4

ENDJOB